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Abstract 



The relation between the notion of crystalline symmetry and characteristic 
time intervals when this symmetry could be observed is analyzed. Several time 



scales are shown to exist for a system of interacting particles. It is only when 
the observation time is much larger than the mesoscopic fluctuation time, the 
notion of crystalline symmetry becomes physically meaningful. The ideas are 



concretized by a two-phase lattice model. 
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X ' 1 Characteristic Time Scales 

Crystalline symmetry is a well-defined geometric notion. But a crystal as such is a 
physical object consisting of many moving particles. In what sense these particles 
form a crystalline lattice with a given symmetry? An intuitive answer to this ques- 
tion would be that it is the average positions of particles which form a crystalline 
lattice, with the averaging performed over sufficiently large time. But, to be math- 
ematically correct, it is necessary to concretize what does mean "sufficiently large" 
with respect to the observation time. The aim of this report is to demonstrate that 
the minimal time scale providing the lower threshold for the observation time is the 
mesoscopic fluctuation time. 

There are several characteristic time scales for a system of interacting particles 
[1]. The smallest typical time is the interaction time 

Tint ~ — , (1) 

V 

in which a s is an effective interaction radius, or effective particle size, or scattering 
length, and v is a characteristic particle velocity, say, sound velocity. When observing 

1 



a system of particles during an observation time T b s that is smaller than Tint, the 
motion of these particles looks absolutely random, with no preferable positions in 
space. 

Another important scale is the local- equilibrium time 

Tloc ~ - , (2) 

V 

where A is a mean free path. The latter is connected with the scattering length a s 
and the average density of particles, p, as A ~ (pa^)^ 1 . The average density p ~ a^" 3 , 
with ao being a mean interparticle distance. Hence, A ~ o^/o? s , from where 

T int ■ (3) 

a s J 

During the local-equilibrium time, particles interact with each other, so that the 
state of local equilibrium is being evolved. 

At a very short time t, one has the dynamic stage, when 

< t < T int , (4) 

and particles move randomly. No spatial symmetry in the positions of particles can 
be noticed for the observation time in the diapason (4). The second time interval 
corresponds to the kinetic stage, when 

Ti nt < t < n oc . (5) 

At this stage, interparticle correlations begin arising due to interactions, but the 
overall motion looks yet rather disorganized, resembling the Brownian motion. In 
solids, ao is just slightly larger than a s , so that the local-equilibrium time (3) is 
only a little longer than the interaction time, t/ oc ~ Tint- The kinetic stage (5) is 
practically absent in solids. One has the hydrodynamic stage, when 

t 1oc < t < oo . (6) 

It is at this stage when the averaged motion of particles is commonly assumed to 
define a crystalline lattice symmetry, under the appropriate external conditions. 

However, in general, it is not the local-equilibrium time (3) that gives the lower 
threshold for the observation time, in order that the lattice symmetry be well defined. 
The actual lower threshold is much larger than t\ oc . There exists one more time scale 
that can be called the mesoscopic fluctuation time 

Tf~ -Tloc, (7) 

ao 



in which If is a linear size of a quasiequilibrium fluctuation [1]. The fluctuation is 
termed mesoscopic because of its intermediate size, being in the interval a$ <C If <C 
L, where L is a characteristic length of the system. Since If ^> a , one has Tf ^> t\ oc . 
In this way, the hydrodynamic stage (6) consists of two qualitatively different 
parts. One is the mesoscopic stage, when 

T loc <t<T f . (8) 

Here, crystalline symmetry cannot yet be rigorously defined for the whole system, 
although the main part of the latter may already be approximately described by a 
kind of symmetry. At the macroscopic stage, when 

Tf < t < oo , (9) 

the averaged motion of particles can be correctly characterized by assigning a con- 
crete spatial symmetry. It is the mesoscopic fluctuation time (7) that defines the 
lower threshold for the observation time which is necessary for making the notion of 
symmetry physically meaningful. In this regard, sufficiently large observation time 
implies 

Tabs > Tf . (10) 

Only averaging over the observation time satisfying the inequality (10) provides a 
mathematically correct definition of the related lattice symmetry. 

2 Two-Phase Lattice Model 

In order to examplify the existence of mesoscopic fluctuations, let us consider a 
lattice model for which such fluctuations not only do exist but their presence makes 
the system more stable. Lattice models describe the physics of solids where an 
essential concentration of vacancies can arise. These vacancy defects appear in 
crystals with high anharmonicity. Among such objects, we could mention quantum 
crystals [2-6], crystallized plasma [7], crystallized white dwarfs [8], crystals near 
liquid-solid transition [1,9,10] and crystals close to vitrification [1,11]. Mesoscopic 
fluctuations of heterophase nature are common for liquid crystals [12]. 

It is worth noting that mesoscopic heterophase fluctuations exist not only in 
crystals but also in other materials. Thus, they are clearly observed by neutron 
scattering experiments in many magnetic alloys [13,14]. Even in such classical fer- 
romagnets as Fe and Ni, there appears local magnetic order in their paramagnetic 
phase, which is also a sort of mesoscopic fluctuations [15,16]. Similar phenomenon 
occurs in high-temperature superconductors where it is often termed phase separa- 
tion [17-22]. More examples can be found in Ref. [1]. 

Now, let us turn to the lattice model to be analyzed in what follows. Consider 
a system of N particles in volume V. Each particle can occupy one and only one 



lattice site of a given crystalline lattice. Not all lattice sites are occupied, but 
some of them are vacant. The distribution of vacancies through the volume is not 
uniform. Some parts of the system contain more vacancies than other, which results 
in essentially different density of particles in the corresponding regions. The location 
and shapes of these parts with different density are not fixed, but fluctuate in space 
and time, forming a system with mesoscopic density fluctuations [23]. Thus, we 
have a heterophase system with coexisting dense and dilute phases. Such a system 
can serve as a cartoon of a crystal with regions of disorder [24] . 

Let us label the quantities related to the dense phase with the index 1 and those 
corresponding to the dilute phase, by the index 2. For instance, the average number 
of particles pertaining to the dense phase is N\, occupying the volume V\, while the 
average number of particles in the dilute phase is N 2 , inside the volume V 2 . The 
total number of particles and the system volume are conserved: 

N = Nt + N 2 , V = V 1 + V 2 . (11) 

Because of the presence of vacancies, the number of lattice sites, Nq, is larger than 
the total number of particles N. Hence, the density of sites is larger than that of 

particles: 

N N 

po> P , po = -^ , p -y- ( l2 > 

By definition, the dense phase has a larger density of particles than the dilute phase: 

Pi>P2, P"=TT (^ = 1,2). (13) 

Introducing the geometric probabilities of the corresponding phases, 

V v 

w u = — , w 1 +w 2 = l , (14) 

the mean density of particles can be written as 

p = wxpi + W2P2 • (15) 

For what follows, it is convenient to define the dimensionless densities 

ra„ = — , n = — = — - , (16) 

Po Po ^0 

for which Eq. (15) reduces to 

n = wirii + W2n 2 . (17) 

A system with mesoscopic heterophase fluctuations is quasiequilibrium and has 
to be described by a quasiequilibrium Gibbs ensemble [1], with a statistical operator 



p(0 ~ e~®^\ where Q(£) is a quasihamiltonian and the set £ = {^(r)} of manifold 
indicators characterizes a distribution of phases in space. This phase distribution is 
frozen at the mesoscopic stage (8), when no lattice symmetry can be defined. Going 
to the macroscopic stage (9) implies averaging over all phase configurations. This 
averaging can be defined [1] as functional integration over the set £ = {^(r)} of the 
manifold indicator functions. In the process of that integration, the notion of the 
heterophase Hamiltonian H arises, defined by the relation 



in which f3 = T" 1 is inverse temperature, ks = 1- From Eq. (18), one has 

H= -Tin /V Q(?) £>£. (19) 

The i-site of the lattice, pertaining to the part filled by the z/-phase, is character- 
ized by the variable e^, taking the value e^ = 1 if the site is occupied by a particle 
and the value e^ = when the site is empty. For the heterophase Hamiltonian (19), 
we obtain 

H = H x © H 2 , (20) 

which is a sum of the phase-replica Hamiltonians 

N ^ N 

H u = w u J2 ( K i - /-0 ei„ + ~ W IY1 $ *i e ^ e i^ > ( 21 ) 

where K^ is kinetic energy; /z, chemical potential; and $y is interaction potential. 
The phase probabilities w u are defined from the minimization of the thermodynamic 
potential 

y = - 1 In Tr «H® , (22) 

in which, for taking account of the normalization W\ + W2 — 1, it is convenient to 
introduce 

w = w\ W2 = 1 — w . (23) 

Then w is given by the equations 



dy d 2 y 



> iSi > ° • ( 24 ) 



In terms of the occupation operators e i; ,, the densities (13) are 



I No 

Pu = ttJ2<^>> ( 25 ) 

v i=l 



where < . . . > implies the statistical averaging, and the densities (16) become 

■y n 
n " = TF J2 <e i"> ■ ( 26 ) 

Recall that, by definition, v — 1 corresponds to the dense phase while v — 2, to the 
dilute phase, so that p\ > p 2 and ni > n 2 . 

It is possible to pass to quasispin representation by introducing the variable 
<y iv = ±1, defined by the relations 

a iu = 2e iv - 1 , e iv = - (1 + o iv ) . (27) 

Then the Hamiltonian (21) takes the form 

y N i N 

H U = U U + -wl J2 ^i3°iv °> - o W » Yl h iv°ii> > ( 28 ) 

8 i#i 2 i=l 



where 



-^ 7V y No 

1 1 "» 



h iu = /J, - Ki - -w u — J2 ®ij ■ 

Z ivn •/ • 

u i+3 

The form (28) corresponds to the Ising model in an external field. 

To accomplish concrete calculations, let us consider the Kac-type interaction 
potential $y = <J>y(iVo), for which 



lim QijiNo) = , lim 

No—too ' Nq^oo 



i JVo 



iVo 



< ex) . (29) 



For this kind of potentials, the problem as is known [25], can be solved asymptotically 
exactly in the thermodynamic limit, when N — > oo, V — » oo, so that N /V — > 
const. In the case of such long-range potentials, the mean-field decoupling 

< (Tij, (7j w > = 5y + (1 - <$jj) < (7^ >< a ju > (30) 

becomes asymptotically exact. 

3 Numerical Results and Conclusions 

With the decoupling (30), we can calculate all thermodynamic characteristics of the 
model. In terms of the occupation operators e iu , this decoupling is equivalent to the 
equality 

< e iv e jv > = - Sij(< e iu > + < e ju >) + (!- 8^) < e iv >< e ju > . 



Calculations result in a system of transcendental equations which we have analyzed 
numerically. Dimensionless densities ri\ = a, n,2 = b, (a > b), as functions of 
temperature T, in units of $ = ^- X^$ij, are presented in Fig. 1, for different 
mean densities. Note that from the condition of heterophase stability, that is the 
second of Eqs. (24), it follows that $ > 0. Fig. 2 shows the probability of the 
dense phase vs. temperature. The shifted chemical potential //*, defined by the 
equations /x* = (pt — K)/Q, K = jj- J2i=i Ki, is given in Fig. 3. The pressure, 
P, specific heat, Cy, and isothermal compressibility, Kt, are shown in Figs. 4 to 6, 
respectively. All curves in Figs. 1 to 6 end at the points where the probability of the 
dense phase becomes zero, after which the heterophase system looses its stability, 
since the phase probabilities, by definition, must be always positive. Finally, in 
Fig. 7, the thermodynamic potential uj w = yT/$ is drawn for the case of the 
heterophase system and also this potential, a>i, for the pure dense system when 
w — 1. As is seen, one always has u w < oo±, which means that the heterophase 
system is thermo dynamically more stable than the pure system. 

In this way, a heterophase system, exhibiting mesoscopic fluctuations, is ther- 
modynamically more stable than a pure system without these fluctuations. At the 
same time, the existence of mesoscopic fluctuations implies the occurrence of local 
instability [1]. In the present case, such instabilities are related to the fluctuational 
grouping of vacancies in some parts of the system, which forms fluctuating disor- 
dered regions locally destroying crystalline symmetry. It is only for large observation 
times, satisfying the inequality (10), when the system displays, on average, a perfect 
symmetry of the lattice, while for times shorter than the lifetimes of mesoscopic fluc- 
tuations, an ideal lattice symmetry cannot be observed. However, such a temporal 
breaking of symmetry makes the system more stable in the long run. 

The results examplified here by a lattice model can be generalized to other phys- 
ical systems. To our mind, the same conclusions concern not only simple physical 
systems but, even to a greater extent, all complex systems, such as societies. Thus, 
a society allowing the occurrence of mesoscopic fluctuations can be called demo- 
cratic, while that one prohibiting any deviation from the prescribed symmetry is 
dictatorial. In the long run, democracy is more stable than dictatorship. 
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Figure Captions 

Fig. 1. The density of the dense phase (solid curve) and of the dilute phase 
(dashed curve) as functions of temperature. 

Fig. 2. The probability of the dense phase vs. temperature. 

Fig. 3. The shifted chemical potential vs. temperature. 

Fig. 4. The dimensionless pressure vs. temperature. 

Fig. 5. Specific heat vs. temperature. 

Fig. 6. Isothermal compressibility vs. temperature. 

Fig. 7. Dimensionless thermodynamic potentials for the heterophase system 
(two lower curves) and for the pure dense system (two upper curves). 
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